Magnetic Molecule on a Microcantilever: Quantum Magneto-Mechanical Oscillations 
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We study the quantum dynamics of a system consisting of a magnetic molecule placed on a mi- 
crocantilever. The amplitude and frequencies of the coupled magneto-mechanical oscillations are 
computed. Parameter- free theory shows that the existing experimental techniques permit observa- 
tion of the driven coupled oscillations of the spin and the cantilever, as well as of the splitting of 
the mechanical modes of the cantilever caused by spin tunneling. 

PACS numbers: 85.85.+j, 75.50.Xx, 75.45.+j, 75.80.+q 



Magnetic molecules exhibit quantum tunneling be- 
tween different orientations of the spin in macroscopic 
magnetization measurements Q . Detection of coher- 
ent quantum spin oscillations, similar to those observed 
in a SQUID [3], would be of great interest. In a crystal 
of magnetic molecules this effect is difficult to observe 
because of the inhomogeneous broadening of spin levels 
and decoherence arising from various interactions. Mag- 
netic measurements of individual molecules would have 
been more promising but insufficient sensitivity of ex- 
isting magnetometers has prohibited such studies so far. 
The effort has been made to observe spin tunneling ef- 
fects in the electron transport through a single magnetic 
molecule bridged between metallic electrodes [i^. 

In this Letter we propose a different approach to the 
detection of quantum oscillations of the spin of a single 
magnetic molecule. It is based upon the resonant cou- 
pling of the spin oscillations to the mechanical modes of 
a microcantilever. The geometry of the proposed exper- 
iment is shown in Fig. [H A magnetic molecule of spin 



magnetic 
molecule 




FIG. 1: Geometry considered in the paper. 



where | ± S) denote two opposite spin orientations along 
the X-axis. A weak ac magnetic field of frequency 
uo = A/fi, applied along the X-axis, will force the spin of 
the molecule to oscillate between the two orientations. 
Conservation of the total angular momentum requires 
that the oscillations of the spin are accompanied by the 
mechanical oscillations of the cantilever (Einstein - de 
Haas effect [El, Q). Consequently, if a; = A/fi coincides 
with a resonant mode of the cantilever, one should expect 
the effect of the ac field on the cantilever. 

If the cantilever rotates by a small angle 54> at the 
location of the spin, the Hamiltonian of the molecule Hs 
becomes 0, [s| 



H's = 



RHsR~^ = Hs+i'^\Hs,S ■5<p = hS-5<p. 



(2) 



Due to strong magnetic anisotropy the molecule can be 
considered as a two-state system. The Hamiltonian of 
such a system, H2, is a projection of the Hamiltonian 
([2]) onto the two states given by Eq. ([T]). These states 
can be viewed as the eigenstates of the Pauli matrix a^. 
For the geometry shown in Fig. [T] the projection can be 
performed by writing 



W2= ^ {^i\{Hs 



with 



a, = |M/+)(M/-|- 

<7, = |* + >(*+|- 



|*_)(*_|. 



(3) 

(4) 
(5) 
(6) 



This gives 



S is deposited on a microcantilever of length L, Fig. 1. 
The y = end of the cantilever is fixed while the y = L 
end is free. The molecule is assumed to be imbedded in 
or firmly attached to the cantilever at = ?/o, with the 
magnetic anisotropy axis being parallel the X-direction. 
Let A be the tunnel splitting of the two lowest energy 
states of the molecule. 



S)) 



(1) 



W2 = -2^<^z 



(7) 



The states with a definite X-projection of the spin are 
eigenstates of da,: | ±5') = -i=(|^+) ± |^-)). The expec- 
tation value of a satisfies the Landau-Lifshitz equation: 



-a X b, 



'eff : 



with 



beff = -Ae^ + 25'A(50^e 



y • 



(8) 



(9) 
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If the cantilever was held stationary {5(l)x = 0), the solu- 
tion of Eq. ([8]) would describe pure quantum oscillations 
of the spin: cr^ = const, ax cx cos{At / h) , (jy ex sm{At/h). 

We are interested in the coupled oscillations of the spin 
and the cantilever. Since the latter contains macroscopic 
number of atoms, its oscillations can be studied within 
a continuous elastic theory that deals with the displace- 
ment field u(r, t) and the local rotation Scf) given by ^] 



(r)=-Vxu(r). 



(10) 



Within such a model the spin of the molecule at a point 
ro can be replaced by the spin field 



i:{r,t)=Scr{t)d{r-ro), 



(11) 



where a{t) satisfies Eq. ([8]). With account of Eq. ([2]) the 

energy of the cantilever becomes 



'He = He 



(fr ht,-{V xu), 



(12) 



where He '^s the part of the elastic energy that is inde- 
pendent of the spin. 

The dynamical equation for the displacement field is 



P- 



dX(3 



(13) 



where aaf3 = is the stress tensor, eaf3 = 

dua/dx(3 is the strain tensor, and p is the mass density 
of the material. This gives 



da 



(E) 

af3 



dXf3 



V X E 



(14) 



where a^^^ = SHE/^eafd- It is easy to see that in the ab- 



sence of the external torque, K^-' = § dAs 



applied to the surface of the body A, Eq. ([T4j) provides 
conservation of the total angular momentum. 



J £r [hT. + p (r X u)] , d3/dt = . 



(15) 



Writing the left-hand side of Eq. ([T4j) in the conven- 
tional form for a cantilever [9] one obtains the elastic 
equation that couples vertical displacements of the can- 
tilever, Uz{y,t), with the oscillations of the spin: 



h^E d'^u^ hS d d 



df^ 12(1-0-2) %4 2Vdydt 



[a^{t)L5{y - yo)] • 
(16) 

Here h and V are the thickness and the volume of the 
cantilever, respectively, E is the Young's modulus, and a 
is the Poisson coefficient, —1 < a < 1/2. 

It is convenient to switch in Eq. (p!6|) to dimensionless 
variables Uz = Uz/L,y = y/L,t = tu, where 



' Eh^ 
12p(l-a2)L4 



(17) 



determines the scale of the eigenfrequencies of the os- 
cillations of the cantilever. By order of magnitude v ~ 
VshjL? where Vq ~ yE/p is the speed of sound. In terms 
of these variables Eq. ([16]) becomes 



d'^Uz d^Uz 



where < < 1 and 



e = 



hS 
ML'^v 



hS 
M 



12p{l - cr2) 



Eh^ 



(18) 



(19) 



is a dimensionless small parameter. By order of magni- 
tude, e ~ hS/{Mvsh), where M = pV is the mass of the 
cantilever. For, e.g., a molecule of spin 6* = 10 on a can- 
tilever of dimensions lOOnmx lOnmx Inm the parameter 
e should be of order 10~^. Eq. ([18]) has to be solved with 
the following boundary conditions: 



duz 

Uz = , -—zr =0 at ^ = , 



d^Uz 
dy'^ 



0, 



d^Uz 
dy^ 



at y = l. 



(20) 



The first two conditions correspond to the absence of 
displacement and the absence of bending of the cantilever 
at the fixed end, while the last two conditions correspond 
to the absence of torque and force, respectively, at the 
free end [9]. 

For the free oscillations of the cantilever (e = 0) one 
writes 

Uz{y, t) = u{y) cos{Ot) . (21) 
Substitution into Eq. ([T8|) with e = then gives 



dy^ 



Hi U = , 



(22) 



Solutions are 



u{y) = {coshz + coshA>:) [cos{hzy) — cosh.{hzy)] 
-\-{smhz — smhhz) [sm{tzy) — smh.{hzy)] . (23) 

The third of the boundary conditions ([2Q|) provides the 
equation. 



cos K: cosh + 1 = , 



(24) 



for the frequencies of the normal modes of the cantilever, 
On = i^n (measured in the units of v of Eq. ([T7|) ). The 
fundamental (minimal) frequency is coi ^ 3.516. The 
next two frequencies are Cd2 ~ 22.03 and Cos ^ 61.70. 
The profiles of the oscillations of the cantilever for three 
normal modes (n = 1, 2, 3) are shown in Fig. O 

To consider coupled oscillations of the cantilever and 
the spin of the molecule we first neglect dissipation and 
write for the displacement 



(25) 
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n=l 



co.t = n 



OJit = 



O.O 0.2 0.4- 0.6 O.S Vl-O 




0.0 0.2 0.4 0.6 0.8 yi.o 




ujn^ respectively, unless these two frequencies are very 
close to each other. The latter can be achieved by, e.g., 
changing A with the help of the dc magnetic field per- 
pendicular to the anisotropy axis of the molecule. At 
A = hwn one should observe the splitting of the mechan- 



ical mode of the cantilever, 



into two modes 



2 



S = 



2Un 



(31) 



The remarkable property of Eq. (|3T]) is that it has no 
free parameters. For a chosen resonance A = hjn^ the 
relative splitting 5 depends only on the position of the 
molecule on the cantilever. This dependence is plotted 
in Fig. El 



FIG. 2: Profiles of the oscillating cantilever at different mo- 
ments of time for n — 1, 2, 3. 



where Rm{t) are functions of time to be determined 
and Um{y) is a normalized eigenfunction ([23|) that corre- 
sponds to the eigenvalue given by Eq. ([24j) . 



/ dy Um{y)Un{y) = Smn • 

Jo 

Substitution of Eq. ([25]) into Eq. ^ gives 



(26) 



E 



dp 



-UJrn^m Um{y) 



e d d 
2 di dy 



[^x{t)5{y-yo)] , 
(27) 

where we have used Eq. ([22|) . Multiplying both parts of 
this equation by Un{y) and integrating over y from to 1 
with account of Eq. ([26]). one obtains linear second-order 
differential equation for Rn{t), 



d^R, 



da. 



where <(^o) = {dun/dy)y=y^. 

Coupled magneto-mechanical oscillations near the 
ground state correspond to small (Jy^ R^ and cjz ^ 1- 
This requires temperatures ksT <^ A. In this case Eq. 
dH]) gives for ax 



"^'^^ '•AV. = 5A<(^-o)-^^^ 



dp 



dt 



(29) 



where A = A/{hiy). Substituting into Eqs. ^ and ^ 
ax{t)^ Rn{t) oc exp(icjF), one obtains the following equa- 
tion for the eigenmodes of the coupled oscillations: 

[co^ - coDiCo' - A') = ^eSAu'^{yo)Co' ■ (30) 

Due to the smallness of e, oscillations of the spin and the 
cantilever occur independently at frequencies A/h and 
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FIG. 3: Color online: The dependence of the relative split- 
ting of the first three cantilever modes on the position of the 
molecule. 



We shall now demonstrate that the above experiment 
can be performed by studying the response of the system 
to a weak ac magnetic field, B(t) = BoGx sin(cc;t), applied 
along the anisotropy axis of the magnetic molecule. In 
the presence of such a field Eq. (|9]) becomes 

beff = - Ae^ + 2SAS(l)xey + bex sin(cjt) , (32) 

where b = gfisBo, with g being the gyromagnetic factor 
for the spin and fiB being the Bohr magneton. To obtain 
the amplitude of the forced oscillations of the cantilever 
we need to include dissipation in the equations of motion 
of the spin and the cantilever. The damping modifies the 
Landau-Lifshitz equation [10] : 

h& -a X beff + 2Q~^a x (cr x beff) . (33) 

Here Qs is the quality factor of the spin oscillations. 
Small oscillations (|cra.^^| <C l,^^ ^ 1) now satisfy 



da± 
"dT 



A {Q-^ ±i)a± = SAu'^{yo)Rn T ib sm{cdt) (34) 



where a± = ax ^ iay and b = b/ {hu). 
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Dissipation of the mechanical motion of the cantilever 
can be introduced by adding the first time derivative of 
Rn to Eq. ([28]): 

(PRn On dRji . -2 r) ^ -/ ^ - ^ d , . , . 

= ^^^^ 

Here Qn is the quality factor of the oscillations of the can- 
tilever at the eigenfrequency uOn- One can now obtain the 
time dependence of Rn and (j± by solving together Eq. 
(|34|) and Eq. ([35|) . The displacement of the cantilever at 
a point y is given by Uzit^ y) = Rn{t)un{y) • The simplest 
way to get the solution is to replace i sm{Cdt) in Eq. (|34j) 
with exp{iCdt) and solve the resulting three linear alge- 
braic equations for Rn^cr± ex exp(icDt). The applicability 
of the formulas obtained that way is limited by the range 
of parameters that provides the condition \(j± \ <C 1 used 
to derive Eq. from Eq. ([8|). It is easy to see from 
Eq. that this requirement is violated for (j_ when 
CO is close to A: The strong pumping of the spin excita- 
tions by the ac field at a; ^ A leads to the breakdown 
of the linear approximation for the dynamics of the spin 
governed by Eq. (|8|). 

At = cjn 7^ A, in the practical range of the quality 
factors, 1 <C Qn,s <^ one obtains 



\uz{yo,y)\ = 



beQn/^\u'^{yo)\\un{y)\ 

20n\0l-A^\ 



\^±\ 



\On±A\ 

(36) 

The parameter A can be controlled by a dc magnetic 
field applied perpendicular to the anisotropy axis of the 
molecule. The condition |cr_| <C 1 determines how close 
to the double resonance, A = cjn, one can use Eq. ([36|) : 
\(jn — A| ^ 6/2. Substitution of \On — A\ ^ b into the 
first of Eqs. (|36|) provides the estimate for the maximum 
amplitude of the oscillations of the cantilever at A ^ cj^- 
max |ix^(^o, ^)| ^ ^\Un{yo)\\Un{y)\ ' (37) 

It is a function of the distance y from the fixed end of the 
cantilever, parameterized by the position of the molecule 

yo- 

For, e.g., n = 1 the parameter S of Eq. ([3T]) should 
be generally of order A cantilever of dimensions 

lOOnmx lOnmx Inm, carrying a magnetic molecule of 



spin 10, will have v ~ 10 s 



10- 



10" 



see Eqs. ([Ill), (ffSl), and jUj). This would give /i = 
3.52i^/(27r) - 30MHz and the splitting, Sfi, of the first 
harmonic of the cantilever in the ballpark of a few kilo- 
hertz. The condition for the detection of the splitting is 
Qn,s ^ l/S. For a single magnetic molecule on a micro- 
cantilever the discrete character of the cantilever phonon 
modes and the absence of the inhomogeneous broaden- 
ing of the spin mode (usually present in a bulk molecular 
magnet) should provide a high spin quality factor. Very 
high quality factors (up to 10^) have also been reported 



for microcantilevers 11, 12[. When the molecule is near 
the tip of the cantilever, Eq. (|37|) with e ~ 10~^ and 
Qn 10^ gives for the tip: \uz\ - 0.005. For a 100-nm 
long cantilever this would correspond to the oscillations 
by half a nanometer. Such a displacement can be de- 
tected by tunneling or force microscopy, as well as by op- 
tical and electrical methods. Working with even smaller 
cantilevers would provide higher values of e and would 
allow to relax the requirement on the quality factor. 

In Conclusion, we have shown that driven quantum 
oscillations of the spin of a magnetic molecule can be 
observed by placing the molecule on a microcantilever. 
Since such cantilevers consist of hundreds of thousands 
of atoms, this would be a remarkable example of a macro- 
scopic quantum effect. Our theory has no free parame- 
ters and, therefore, it must be helpful in designing the 
proposed experiment. 
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